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Abstract. The average charged hadron muhiphcities induced by top quark decays are calculated in pQCD 
at LHC energies. Different modes of top production are considered. Proposed measurements can be used 
as an additional test of pQCD calculations independent on a fragmentation model. 

PACS. 14. 65. Ha Top quarks - 12.38.Bx Perturbative calculations - 13.85.Hd Inelastic scattering: many- 
particle final states - 13.85.Ni Inclusive production with identified hadrons 



1 Introduction 

The study of unstable heavy particles like W, Z bosons, 
top quarks and others (arising in different extensions of 
the Standard Model) is one of the leading directions in the 
modern high energy physics. To determine particle param- 
eters (charge, mass, width, decay modes etc.) we have to 
look deep inside their production and decay mechanisms. 

In this article processes with top production are con- 
sidered due to its specific properties. Top quark is ex- 
tremely elusive object. Because its mass is so large (m^^^ = 
172 ± 0.9 ± 1.3 GeV) [T], it can decay into on-shell W- 
bosons, i.e., the two- particle decay mode t — >■ bW^ is 
kinematically possible. The SM predicts the top quark 
to decay almost exclusively into this mode. The on-shell 
W-boson can then decay leptonically or hadronically with 
coupling strengths given by the Cabibo-Kobayashi-Maskawa 
(CKM) matrix. The top quark decay proceeds extremely 
fast, in less than tj = l/ri ~ 5 x 10~^^ s, which is shorter 
than the time scale to form hadrons T^ad - 1/Aqcd « 
3 X 10~^^ s and almost 13 orders of magnitude smaller 
than the lifetime of hadrons, which involve the next heav- 
iest quark, Tb — 1.5 x 10~^^ s. The width Ft acts as 
a physical "smearing", and the top production becomes 
a quantitative prediction of pQCD, largely independent 
on nonpcrturbativc phcnomcnological algorithms. That is 
why top measurements are directly related to pQCD tests. 

The principal observables used in tests of pQCD are 
typically measurements of jets, high transverse momen- 
tum particles and event shapes. There arc relative advan- 
tages and disadvantages in using these observables. Jet 
measurements are expected to have a close correlation in 
direction and momentum with the parton which gave rise 
to it. However, several complications concerning jet def- 
initions (algorithms) arise when using jet measurements. 
While the infrared and collinear safety of event shapes al- 
lows safe perturbative predictions, we have to take into ac- 
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count resummation of large logarithms. Evolution of struc- 
ture functions can be predicted by pQCD calculations, but 
we have to consider different kinematical regimes (see, for 
example, recent review on QCD tests 0). The transition 
from partons to hadrons cannot be accomodated within 
perturbative QCD. Fragmentation mechanism is usually 
simulated by the use of additional approaches (string frag- 
mentation, cluster fragmentation) . It was shown that 
measurements of average charged hadron multiplicitiy in 
a jet (especially jet produced in heavy quark decay) can 
serve as a precise test of pQCD independent on a frag- 
mentation model. 

Measurements of quark and gluon jets show visible dif- 
ferences between them. For example, gluon jets are fatter, 
softer and have higher multiplicity [S]. In events induced 
by heavy quark jets the hadron multiplicity is smaller than 
in analogous events triggered by light quark jets [B]-[5]. 
The situation is similar to the classical theory, where the 
more is the mass of a charged particle, the less intensive is 
the radiation from it. The dependence of the multiplicity 
in a quark jet on the " primary quark" mass was observed 
in previous experiments, in particular, at the SLC, KEK, 
LEPI and LEPII (see [S] and references therein). Leading 
order pQCD predicts "specific scaling" I3],[S]-[II], i-e. en- 
ergy independent difference between average charged mul- 
tiplicities of light and heavy quark jets. The advantage 
of multiplicity measurement is that we do not need large 
number of events, and it can be used in rare processes like, 
for example, single top production. 

The QCD radiation associated with tt production has 
been treated earlier in [T^ll3j . These papers consider the 
effect, that gluons are also radiated from the b's from t- 
decay. This effect has been also taken into account here. 
Recently top production at the ILC e+e~ collider was con- 
sidered [3] . Following the basic idea of Ref . [1] it is possible 
to analyse the case of top production at the LHC. Since in 
pp events the energy of parton interaction is not fixed and 
changes in the range from about 400 GeV up to several 
TeV, we can study the evolution of average multiplicities 
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Fig. 1. Parton level processes of top production at LHC con- 
sidered in this paper: a) dominant tt production in the gluon- 
gluon fusion; b) s-channel single top production; c) t-channel 
single top production; d) tW production. 



in this energy domain. Of course, several obvious difficul- 
ties can arise (color interconnection, interference effects). 
It is shown, that there are some cases where we can suc- 
cessfully avoid the above complications. 

The paper is organized as follows. In the first section 
there is an outlook of the basic model. The second one 
is devoted to numerical results for different mechanisms 
of the top production at the LHC. Complicated formulae 
and calculations are collected in Appendices. 



2 Model 

The model used in this section is analogous to the one 
presented in [3]. In this article we consider processes of 
the type p + p^tt-^Xorp + p^t + X (see Fig. [T]). 
The basic formula for inclusive production of the system 
M in the collinear approximation looks as follows 



dxidx2d<P 



M 



= y2Mxi).fj{x2) 



da 



d^ 



M 



{xiX2S\ {^A/}), 
(1) 

where fi{x) is the probability to find parton i (quark, anti- 
quark or gluon) with the longitudinal momentum x^fsjl 
in a proton. Renormalization and factorization scales are 
hidden in / and (j. Here M is the system of different fi- 
nal states like ti, t6, ig, tW , which corresponds to dif- 
ferent mechanisms of the inclusive top production at the 
LHC. X includes beam remnants (A'beam) and also the 
secondary radiation induced by color interactions inside 
M [Xm) plus possible interaction between M and beam 
remnants {XM-beam), if M is not a color singlet. Vari- 
ables a; 1.2 are fixed in every separate event and can be 
calculated experimentally. Usually the sum in ([T]) can be 
approximated by a single factorized term which includes 
parton distributions multiplied by the amplitude squared 



of parton-parton cross-sections: 

g + g if + (dominates at the LHC), (2) 

q + 4 ^W* ^ib + XM, (3) 

h + q^U + Xm, (4) 

6 + + (5) 

q + g^tbq' +Xm, (6) 

g + g^tbW + XM, (7) 

where q, q' denotes corresponding light quarks. For our 
purposes it is enough to consider initial parton collisions 
instead of pp process, since we have to calculate only 
the charged hadron multiplicity of the M plus Xm in 
a separate event. Complications concerning Xteam and 
XM-beam are discusscd below. Xm comes from virtual 
gluon radiation. 

For the average multiplicity of hadrons in M -I- we 
can use the expression similar to Eq. (5) in [3]: 



Nti+x,AQ') = UM + 
d'^k 



ab I 



(2^)4 ",.(91, 92, k) d^^:,ik) dl^,{k) n'il (fc) , (8) 

where d^l^{k) = i d^^{k)/{k^ + iO) is the propagator of 
the gluon with momentum fc, qi^2 are momenta of initial 
partons. Here and below (a, b) and (a', b') denote color in- 
dices, Q = \/ (^1 + 92)^ is the energy of colliding partons. 

The first term in the r.h.s. of Eq. (|5]), 7ij\/, is the multi- 
plicity from the fragmentation of leading particles in the fi- 
nal state. For example, in the process ([2]) um = tlu = 2nt, 
where nt was calculated in [3]. In other processes is 
appropriate combination of multiplicities which are taken 
from the analysis of data and pQCD calculations: 



.a'b' , 



nt{t -> hadrons) = 41.03 ± 0.540], 



(9) 



nt{t Ivi + hadrons) = 21.9 ± 0.53[3], (10) 

(11) 
(12) 



nw{W hadrons) = 19.34 ± O.lOg 
= 2.6, Ub = 5.5[IU], = 1.2[T3]. 



Tensor 



'^p [k) is given in |3]: 



(13) 



where dimensionless quantity ng{k'^) describes the average 
multiplicity of hadrons in the gluon jet with the virtuality 
k^. It is, of course, gauge invariant, and depends only on 
the virtuality fc^. 

The quantity ng{k'^) cannot be calculated perturba- 
tively. It is usually assumed that the average hadron mul- 
tiplicity is proportional to ng{k^, Qq), i.e. the average mul- 
tiplicity of (off-shell) partons with the "mass" Qq (the so- 
called local parton- hadron duality): 



ng{k^) = ng{k',Ql)K{Ql) , 



(14) 



where K{Qq) is a phenomenological energy-independent 
factor. The QCD evolution equations for both ng{k^,Q'^) 
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and 



dp 



p 



2■"^(^•^<^o) 



(15) 



are derived in [3], and also in old works [TSllTB] . We use 
the "conventional standard" value Qo = 1 GeV for fur- 
ther numerical calculations. Here Ng(k^,Q^) is the av- 
erage multiplicity from the gluon jet whose virtuality 
varies up to fc^. Very often Ng is erroneously called the 
average multiplicity of the gluon jet with fixed virtuality 
. This meaning should be addressed to rig only. 
In this paper we will use two phenomenological expres- 



sions (to estimate theoretical uncertainties) for 
can be found in [3]: 



n'g{e)^e—Nl{k\Ql), 



iVi(fc2) =3.89 + 0.01 exp 



1.63 




/li =0.87 GeV(QCD motivated), 
iV|(fc2) =4.21 + 0.012 In^ 72 ' ^2 = 0.93 GeV 

^2 



which 
(16) 

(17) 
(18) 



Parameters for these functions were obtained in [4] (see 
eq.(3) and Fig. 8 in this reference, where Ug corresponds 
to the function Ng in the present paper) by the fitting of 
the data from [TT] . 

In Eq. ([5]) the first factor of the integrand is given by 



7T"^(gi,(?2,fc) = 



n 



n 



{aihi }; ah 



(19) 



where 7T|°'^'|' ^^(gi, 92, fc) can be calculated in the first 
order in the strong coupling constant as the amplitude 
squared of the corresponding process ©-([T]) with Xm = g 
normalized to the total rate of the process without Xm- 
Sum in the product of polarization vectors of initial gluons 
{Nin.gi. = 0, 1, 2) forms usual factor (here we use the axial 
gauge since it simplifies much theoretical calculations) 



Pi 



d>'^'(q.,n)= J2 
A=l,2 



Pi O 

7li % 



qn qn^ 

0, = 0, n is an appropriate four- vector 
in the corresponding process (see Appendices B,C). 



(21) 



where e^^).^g,;p 



The quantity ([T^ satisfies the equality 
k^^nf,{q,,q2,k)^Q 



due to the general theorem p/7j: If the QCD amplitude 
is written 

Aqcd = 



then one gets zero if any number, > 1, of the polariza- 
tion vectors £pj{kj) and/or ej^. (k^) are replaced by kj^^^ 
and/or Ki ^ i,. respectively, provided that all these kj 's and 
Ki 's, with the exception of at most one of them , satisfy 
k^ = and = 0. 

If we take into account Eq. (PT|) and introduce the 
function 

n{Q^e,kq,,kq2) = {-gn5abnf,{Q\k\kq,,kq2), 

(23) 

the final formula for the multiplicity will look as follows: 

d'^k 



nu + 



(27r) 



-n{Q\k\kq^,kq2) J^Ngip^) 



(24) 



Concrete form of the function U for different processes 
can be found in Appendices B,C. 



3 Numerical results of calculations 

In this section we consider numerical results for average 
charged multiplicities in different processes of top produc- 
tion at the LHC. Below we consider the phase space when 
final jets have low transverse momentum cuts Pt , and the 
final gluon jet can not be experimentally separated from 
one of final quark jets (i.e. gluon jet lies within the cone 
cosOgq > R = 0.9, where 9gq < 0.45 is the angle between 
the gluon and quark jets). 

Let us begin with the inclusive 1 1 production ^ . The 
total cross-section of the inclusive process pp ^ tt + X is 
about 833 pb at 14 TeV. In this article we consider only 
the gluon-gluon fusion mechanism of this process since at 
LHC it is dominant. Numerical values for N^^ are given 
in the table [T] and on the Fig. [2] Here and below the- 
oretical errors are estimated by the use of two different 
parametrizations (|17[) .(|18 p for the hadronic multiplicity in 
a gluon jet. The average charged multiplicity in different 
decay modes (hadronic, semileptonic and leptonic) can be 
calculated as follows 



^ti^hadronsiQ) ' 
tt^Wi+hadronsiQ) ' 
^tt^l+l- UiPi+hadronsiQ) ' 



:2n 



H 



--2ni 



+iVf,(Q), 

mQ). 



(25) 
(26) 
(27) 



As you sec on the Fig. [21 the dependence of N^^ on the 
energy is visible. In this work we assume that color re- 
connection of tt and beam remnants is small due to the 



Table 1. Multiplicity A^^^^ for different cuts of jet transverse 
momenta Pt and the energy of gluon-gluon collision. 





Q, GeV 


Pt, GeV 


600 


1000 


1500 


10 


2.82+0.07 


8.29+0.2 


15.75+0.36 


30 


0.76+0.02 


2.96+0.06 


6.58+0.12 


50 


0.4+0.01 


1.3+0.02 


3.63+0.05 
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Fig. 2. tt production. Multiplicity N^^ versus Q for different 
cuts of jet transverse momenta. Top-down: Pt — 10 GeV — > 
Pt = 30 GeV ^ Pt = 50 GeV. 
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Fig. 3. S-channel single top production. Multiplicity N^j^{Q) 
versus Q for different cuts of jet transverse momenta. Top- 
down: Pt = 10 GeV Pt = 30 GeV ^ Pt = 50 GeV. 



strong suppression of this processes with high transverse 
momentum transfer (typical jet transverse momentum cut 
at the LHC is about 20-40 GeV). Also t and t fragment 
independently after the interaction inside the tt system. It 
looks similar to the process of W~^W^ fragmentation in 
e"'"e~ annihilation. The effect of possible color reconnec- 
tion was investigated by comparing hadronic multiplici- 
ties in e+e^ W+W~ qq'qq' and e+e" ^ W+W^ ^■ 
qq'lvi. No evidence for final state interactions was found by 
measuring the difference < n'l^ > — 2 < nl^qi^, > [TB].p^. 
The values for average charged multiplicities can be com- 
pared with the present LHC data on the inclusive tt pro- 
duction. 

The case of s-channel single top production ([3]) is close 
to the e~^e~ one, since the final state is a result of W decay, 
i.e. color singlet. That is why we have no color reconnec- 
tion with beam remnants. However, the cross-section of 
this process is rather small (about 11 pb at 14 TeV), and 
the experimental task on the extraction of the multiplicity 
looks more difficult than, for example, in t-channel single 
top or tt production. Numerical values for N^j^ are given 
in the table [2] and on the Fig. [3] The average charged 
multiplicity in different decay modes can be calculated as 
follows 



The process of t-channel single top production pp — > 
t + X has higher rate (about 245 pb at 14 TeV) than the 
previous one, but we have to make the same assumptions 
concerning fragmentation and color reconnection processes 
as in tt production. Here calculations for the parton level 
process (jH) are presented. Numerical values for Nf^, are 
given in the table[3]and on the Fig. 01 The average charged 
multiplicity in different decay modes looks as follows 

K'^hadronsiQ)=nf + + 7Vf,,(Q), (30) 
K'^l^.+hadrons{Q)^nf + n, + Nf^,{Q). (31) 

As you can see on the Fig. |4j the value of , is rather 



Table 3. T-channel single top production. Multiplicity N'^^, 
for different cuts of jet transverse momenta Pt and the energy 
of parton-parton collision. 





Q, GeV 


Pt, GeV 


600 


1000 


1500 


10 


6.23±0.18 


7.65±0.22 


8.62±0.24 


30 


2.4±0.07 


2.77±0.075 


3.29±0.08 


50 


1.32±0.038 


1.59±0.04 


1.76±0.044 



N. 



tb->hadro7is(Q) — '"'t + '^fc 



(28) 

(29) 



10 r 



The energy dependence is not so strong as in the previous 
case (see Fig. [3]). 



Table 2. S-channel single top production. Multiplicity N^^^ 
for different cuts of jet transverse momenta Pt and the energy 
of parton-parton collision. 





Q, GeV 


Pt, GeV 


600 


1000 


1500 


10 


11±0.32 


14.8±0.41 


18.7±0.5 


30 


6.55±0.19 


10±0.27 


13.2±0.33 


50 


4.33±0.12 


7.55±0.2 


10.4±0.25 



600 



500 



1000 



1200 



1400 



Q, GeV 



T-channel single top production. Multiplicity Nj^^,{Q) 



Fig. 4 

versus Q for different cuts of jet transverse momenta, 
down: Pt = 10 GeV ^ Pt = 30 GeV ^ Pt = 50 GeV. 



Top- 
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Table 4. Multiplicity N^^y for different cuts of jet transverse 
momenta Pt and the energy of parton-parton collision. 





Q, GcV 


Pt, GeV 


600 


1000 


1500 


10 


2.54±0.06 


5.83±0.14 


10.27±0.23 


30 


0.85±0.017 


2.27±0.04 


4.49±0.076 


50 


0.42±0.008 


1.25±0.019 


2.6±0.036 




Q, GeV 



600 800 1000 1200 1400 



Fig. 5. Multiplicity Nf-^iQ) versus Q for different cuts of 
jet transverse momenta. Top-down: Pt = 10 GeV Pt — 
30 GeV Pt = 50 GeV. 



small in the wide kinematical region, and energy depen- 
dence is not strong. It is important for the estimation 
of the multiphcity from beam remnants plus color 
reconnection effects, since values , , Ug arc fixed 
by previous measurements and Nf^,[Q) <C rit. From this 
point of view the t-channel single top production looks 
the most interesting process for the multiplicity measure- 
ments. 

tW production has intermediate cross-section of the 
order 62 pb at 14 TeV which lies between s- and t-channel 
single top production rates. Probably, specific signature 
of this process would help in the measurements proposed 
in this work. Numerical values for Nfy^ are given in the 
table [4] and on the Fig. [5] The process ([5]) has 3 decay 
modes. The corresponding average charged multiplicities 
are 

N]^W^hadrons{Q)^n^ + + Nf^{Q), (32) 

N^w-.in-.,,,+hadrons{Q)-^n]; + iVf^(Q). (34) 

The energy dependence is also visible and can be used to 
test QCD calculations. 



4 Discussions and conclusions 

In this article we consider four processes with top produc- 
tion at the LHC. Average charged hadronic multiplicities 
were calculated in perturbative QCD. Hadronic multiplic- 
ity in a gluon is fixed by low energy data. 



There are several important tasks that could be solved 
by multiplicity measurements: to test QCD calculations 
independently on fragmentation models, to check in- 
dependent fragmentation of heavy quarks, to check parton- 
parton CM. energy dependence of hadron multiplicities, 
to estimate multiplicity from beam remnants plus from 
color reconnection effects in t-channel single top for fur- 
ther use in other processes. We can calculate also the dif- 
ference ANqq = Nq — Nq to cancel effects of color recon- 
nection and beam remnants. 

There are several assumptions in the present work: 

— independent fragmentation of on-shell top quarks in ti 
production; 

— color reconnection effects in the interaction of jets with 
beam remnants (for nonsinglet production of ti, tq' , 
tW) are suppressed for large lower cuts in jet trans- 
verse momenta. As you can see on the Fig. [51 the 
fragmentation pattern in tt production is rather com- 
plicated. We have beam remnants with low transverse 
momenta interacting with jet remnants with large trans- 
verse momenta. Amplitudes for such processes are sup- 
pressed for high Pt since they are propotional to the 
inverse power of the momentum transfer squared t in 
the parton-parton interaction, and 

i= (pr-^ - p'o^'r - {p'r-- - iff? - 

where p'"^''™ and p^^* are momenta of beam and jet 
partons correspondingly, and from the kinematics 

p^-™ ~ \pT1 » Pi'' ~ bfl, pf* » pT''. 




Fig. 6. Complicated fragmentation pattern of the inclusive 
tt production in pp collisions, a) beam remnants; b) hadrons 
arising from the color interaction of beam remnants with the 
final state radiation from quarks (suppressed for high-Pt of 
final quarks or gluons); c) hadrons from top fragmentation; d) 
result of t*t* interaction. 
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We arc interested only in c and d types of fragmenta- 
tion depicted on the Fig. |6l 
— in the s-ehannel tb singlet produetion there is no color 
reeonneetion with beam remnants. 

All the above assumptions are based on low energy data 
and theoretical estimations of pQCD and also can be 
checked at the LHC. 

From the experimental point of view t-channcl single 
top production (j?]) is the most convenient case, since the 
energy dependence of the average charged hadronie multi- 
plicity is weak. We can estimate quantitatively effect 
of color reeonneetion of jets with beam remnants 
to cheek our assumption on its suppression. Then we can 
use this estimation to improve our predictions for other 
channels of top production. For this task it is also useful 
to extract multiplicities in different decay modes of top 
quarks. 

The final experimental task is to extract number of 
tracks in jets which are produced in top quark decays. 
To estimate experimental efficiencies and dependence on 
a fragmentation model we can use any MC generator for 
top production. At the same time with the top-mass recon- 
struction procedure (in hadronie mode) we could extract 
number of tracks which are included into hadronie cluster 
from single top or top anti-top decays. At the moment we 
have a good chance to make the new independent test of 
QCD by the use of recent LHC data at 7 TeV. Other ex- 
perimental aspects of such measurements will be discussed 
in futher works. 



Appendix A 

Let us consider the process 

partoni{qi) + parton2{q2) ^ c{pi) + d{p2) + g{k) (35) 

and put — Q. k'^ = > (we put also mi, = 
for processes (O,®, since corrections are of the order 
ml/ml <^ 1). In the CM. frame of colliding partons we 
can write: 

gi = |(l,0,0,l), 92 = 1(1,0,0,-1), 



g = 91 +92, ^ = (91 - 92)/2, 



(36) 



, qk VD VD 



I 



2Q 

pi ^ I - k, p2 ^ q - I 



1 1 sin^i COS0, \l I sin^i sin0, |Z | cosfl; 



(37) 



D = {qkf 



K\ \l\ 



/ 4g2 

(9 - k f = Q'^ + K^ - 2qk, Z = f 



m 



2, 



(38) 



cos Oki = cos &!; COS Oi + sm sm Ui cos or 

cos 61 — cos 6ki cos 6k + sin 9ki sin 6k cos cj)* . (39) 



After change of variables phase space looks as follows: 

d^k dH 



dK^ / d^2 



dK' 



(2^)4 (27r)4 

{2tt)5{pI ~ ml){2Tr)S{pl - ml){2Tr)S{K'^ - kf.k'') 
dK^ d{qk) d cos 6k 

dZ d(j) d cos 61 



2(2^)2Q2 

s{n{z)) 



2(\-kY^D 



1 



Ck, 



2(47r) 



3n2 



dK'' 



d{qk) 



dcost 



dZ 



d<j) 



(qk)^ Ck,- Z- 

T-L{Z) = cos 6ki (Z) — cos 6k cos 9i — sin 6k sin 61 cos <f>. (40) 

Here we keep the integration in since the gluon is vir- 
tual. 

Then we have to cut jet transverse momenta from be- 
low to suppress color reeonneetion with beam remnants 

Pi.± > Pt,i, k± > Kt, 



P.,3<^pl-Pt% k,<^k^-Kl 
In this paper Kt = Pt,i = Pt- Finally we have conditions 



|cos 0,11^1 < Jl'-Ph. 



(41) 



||Z|cos6l, - |fc|cos6lfc| < - ky - P^^. (42) 

For limits in the above integrals without conditions (|1T|) . 
we can write 



■\J (mi + m2f 



Kt 



m+ = 



mi A 



mf + Pi 



Ck 



± 



D - Q^Kf 



D 



> ^ [mi + ?7l2)2 + if 2 _^ 



Ql + A'2 > mi +1712 + Qo, 



(43) 



Taking into account the inequality (cos dki can be obtained 
from S{{k — lY — m\)) 



I cos 6ki\ 



Ai 



< 1, 



(44) 
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where 



where 



A, 



and |i|(see ((39|) ) depend on Z, we can obtain limits: 
Z± = 



(45) 



qk 

Di2{s,mi,m2) = 

^ (mi + 7712)^ 



± V^V£'i2(s,mi,m2) 



1 - 



(mi — m2)^ 



(46) 



For multidimensional integration it is convenient to intro- 
duce undimensional variables and make appropriate sym- 
metrization of the function under the integration: 

(j) = X^TT, .T0 G [0, 1], 

s 

CVDy/Di2is,mi,m2), C e [-1, 1], 



cos Ok 



D 



, re [-1,1], 



qk = {qk)_ + 7j{{qk)+ - (qk).), 77 G [0, 1], (47) 

f d(j) 

i d{qk) dcosOk dZ — f{qk, cos 6k, Z, cj}) = 

J TT 



drj dr dC dx^ V f 



fiv, T, -C, x^) + f{i], -r, -C, a;^), 
V^{{qk)+-{qk)^) X 



D - Q^Kfy/Di2is,mi,m2), 



(48) 

(49) 
(50) 



where / is equal to / after the change of variables. 

In this paper we consider the case, when the final gluon 
jet can not be separated experimentally from one of final 
quark jets: 

cos 6*3, > R = 0.9, or Ogq < 0.45. 

The above inequality leads to the following conditions 

cos 9p^k = -cos 9ki > R, (51) 
or cos 9p^k = -cos6'fe/|^.^_^. > R, (52) 

where cos9ki is expressed in terms of variables C, ?7, x^, t. 

Let us denote conditions l|iT|) . . ((5T|) . ([5^ as a prod- 
uct of corresponding 0-functions 



fT,R 



P -P^2-\cos9i\\l\) X 



(l - kf 



p2 



l\ COS 9] — \k\ cos ( 



' (— cos9ki — R) or 6* ( — cos 



R 



(53) 



9{x) 



1, X > 
0, X < 



Now we can rewrite the second term in the r.h.s. of Ref. 
as follows 



(2^ 
d 
dp 



n{Q^,k^,kqi,kq2)fTM X 



32 = fc2 



h 



(54) 



2(4' 
(9fc)+ 



d{qk) J dcOs9k J dZ J ^ fT,Rn2^3, (55) 
Ck.- Z_ 

+c" 



h^- T7. / d cos 61,1 772^2, (56) 



c 



167r 

-c 

12 _ Di2iQ'^,mi,±,m2,±) 
-Di2(Q^,mi,m2) 



(57) 



where 712^.3(712^2) is the amplitude squared of the corre- 
sponding process (121)-® with (without) gluon radiation, 
which is calculated in Appendices B,C. For simplicity we 
put all the coupling constants to unity in these quanti- 
ties. Here Qs = a/IttoJ is the QCD coupling constant. All 
tensors are contracted as in ([T5| . (P^ . 

Different kincmatical invariants of the process psp can 
be expressed in terms of Q^, K^, qk, cos 9k, Z, <f) (and then 

Q^ K^, V, T, c, x^y. 



9192 = — , P1P2 = 



Piqi 
P1Q2 

P2qi 



+ Z - 2qk 



Z-2qk 



4 

Q'^Z 



Al - Ak, 



- AI + Ak, 



Al, P292 = ; + Al, 



Pik 
P2k 

qik : 

Ak = 



2qk - Z - K"^ + mf ~ ml 



'^ + Ak, q2k=^^-Ak, 



■ COS Uk, 



Al = Al cos( 



Bi sin 9k cos c 



(58) 
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Fig. 7. Amplitudes of the process gg — !> QQ. 




Fig. 8. gggQQ amplitudes. 



Q^-qk 



c 



B, 



; " VD ( ml-mn\ 

Di2[s, mi,_L, TO2,_l) — I 1 -z I 



91 

4 



|Z \^-Af = 9^^D,2{s, mi,x, m2,x)v/W^. (59) 



Appendix B 



Here we consider amplitudes for the top-antitop produc- 
tion For the amplitude of the process without addi- 
tional gluon radiation we have three diagrams of Fig. [71 
and I2 can be calculated as follows 



xt ) In 



1 + 



l-'ixt 



3x„, I — 18a::„, — 7 — 16- 



Xm = mt/Q, Xrn± = mt,±/Q. 



(61) 



For the amplitude of the process with additional gluon 
we have five kinds of diagrams (see Fig. [5]): 



7A {qs -pi+mt) d^^iqi + qi) 



ii23 ^ ([213] _ [123;, f 2 . , I ^ ^2 

{ql - 2^193) {qi + q2y 

{{qi - q2)pgp,p2 - (2gi + q2)p2gppi + (2^2 + qi)pigpp2} , 

d^^'ipi +P2) 



([312] + [213] - [321] - [123]) 7A- 



{(2(73 - {Pi +P2))\'9p3P' + (2(pi +P2) - qz)p39\'t3'- 



92 j 



(93 + (pi +P2))/3'5P3A'} . , ^2 

(91 +92)^ 

{(<?! ~ q2)pgpip2 - (2gi + g2)p25/3pi + (2^2 + gi)pi3^ip2} 



123 



7A- 



d^'^(pi +P2) 



(Pl +P2Y 

{([312] + [213]) {gi3p29piP3 + gi3p3gpiP2 - '^gi3pi9p2P3) + 
([321] + [123]) igppigp2P3 + ff,9p33piP2 ^ 25/jp,gp,p3) + 
([132] + [231]) {gpp,gp2P3 + 5/3p25pip3 - 25/3p3ffpiP2)} : (63) 

where wc consider all gluons as initial particles. Then we 
can calculate 



T2^?,°'p\°^P2pSli^l2,qi) 



123 



^132 ^ ^.13 ^ ^.31 ^ ^3 



l213 



231 



l312 



321 



^ (4123 + 4132 +^231) +^123 



1^2 

jj {aibi}:ab 
2—^3 {piCJj};/.ii/ 

-^2^3 P1P2A 



(64) 



;,(gi, 92, -fc)T2%3 ^i^t(gi, 92, -fc). (65) 



Here q^ = —k since li^e /lawe one gluon in the final state 
with momentum k, Af^ = A, ^'^^^^^ (gi, 92, 93,^1,2), 
[ijfc] = t°-H°-H°-^ ^ are SU(3) matrices, it and v are Dirac 
spinors, 7p are Dirac matrices, p = p^^^^. 

If we apply contractions (fTO)) . (j23[) to (|55|) and take into 
account the theorem (it was checked by direct cal- 
culations for (|M)) ) then we obtain g^Il2^3 for the pro- 
cess ([2). Since the final expression for 772^3 is very com- 
plicated, we evaluate it numerically. To get the final result 
for the tt multiplicity induced by gluon radiation (A^^^^ on 

the Fig. [TJa,) we have to substitute I2 and il2-s.3 for this 
process to (|54)) - ([55]) . 



A]^^ ^ u{p2,mt)Al^\>{pi,mt), 
A'^-' [321] 



(62) 



123 _ rQoii^P3 iP2 -qs + mt) 7p2 (gi - Pi + nit) 7pi 
{ql - 2p2q'i){ql - 2piqi) 

A^ = ([321] - [312]) 7P3(P2 -g3 + m07Ari^^(gi+ g2) 

('?|-2p2g3) (gi+92)2 

{{qi - q2)l39ptp2 - (2gi + q2)p29Ppi + (2<?2 + '?l)pi5;3p2} , 



Appendix C 

For simplicity here we set all coupling constants to unity. 
In this section we consider calculations for processes ([3])- 
®. 

Let us introduce some functions for futher calcula- 
tions. One of the functions is the WQq vertex squared 
(see Fig. Hi) 



A'^°^ = u(P2,»712)7q(1 - 75)u(pi,'^i), 



(66) 
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a) 




-£■(1): 2 

' / = Pi aPl a' X 



|a'2 [-Z{2qk - Z) + 2m^{qk - Z) + m'^] 
Z{2qk - Zf - m^iiqk^ -8Z qk + 3Z^) ^ 
m^iiqk - 3Z) - 



~ (1)- 3 

Fig. 9. Amplitudes for calculation of functions T'"-^^ (a) and -^aa' — P2 aP2 a' ^ 



(b). 



SCa{ 9ctot'PlP2 Pi aP2 a' P2 aPl ( 



^^aa.' paPl pP2 a 



(67) 



The next one is the squared amplitude of the process W 
Qqg whieh is shown on the Fig. ^jp. 



7a (1 - 75) ( -Pl - + mi 7^ 



10 + 2pik 



7p [P2 + k + m2 7q(1 - 75) 



v{pi,mi), (68) 



+ 2p2k 

^L'l' ,,,iPi,m,,p2,m2,k,K) = CACpAi'lAiV;, (69) 



(1) 



aa' fj,fj. 



,{pi,m,p2,0,k,K) -cr + 



k^'k^' 



128 



{K^ + 2pikf {K^ + 2p2kf ^ 



(70) 



where 



T-(l); 1 _ 
■^aa' — 9aa' X 



[Z{2qk - Z) + 2m^{Z - qk) - m^] + 

g2 \^2K'^ (ni^ +qk- zf - 2qk{2qk - Z)Z~ 

rr?{Aqk~ Z)Z + 27ri^{qk- Z) + m^ + 

isT^m^ [2qk'^ - 2Z qk + Z^ + 2m^{qk - Z) + m* 
Z{2qk - Z){2qk^ - 2Z qk + Z^) + 
2i-n?{-2qk^ + %qk^Z - IZ'^qk + 2Z^) + 

2m'^{-2,qk^ + 7Z qk - 3Z^) - Am^{qk - Z) ~ m!" \ {71) 



[-Z{2qk - Z) + 2m^{qk - Z) + m^] 
Z'^{2qk ~ Z)- m^Ziiqk - 3Z) + 
m^{2qk-3Z) + m^ 



^aa' = (Pl aP2 a' + P2 aPl a') X 

[~Z{2qk - Z) + 2m^{qk - Z) + m^] - 
[~Z{2qk - Z) + 2rr?{qk - Z) + + 
qk{2qk - Z)Z + 2m^qk^ + m^qk^, 

la' ^ (Pl oka' + ^qPI Q') X 
Q2 

[-Z{2qk -Z) + 2m^{qk - Z) + m^] 



2qk^ 



K'^ [-Z{qk - Z) + m^{qk - 2Z) + m^] H 
^Z{2qk - zf - m^{qk - 2Z){2qk - Z) 



— {6qk^5Z)^m'^ 



(1); 6 



(P2 Q^Q' + kaP2 Q') X 



1^ {-Z(2qk - Z) + 2m^{qk - Z) + m^] 
K'^ [2qk^ -3Z qk + Z^ + 7n^{3qk - 2Z) -\ 
^Z^{2qk - Z)+ m^{Aqk^ - 5Z qk + 2Z^) 



-{Sqk - hZ) + m' 



(1); 7 



aa'apPl P2 ^ 

Q2 [-Z{2qk - Z) + 2m^ {qk - Z) + m^] 

K'^ [2qk^ - 3Z qk + Z^ + m^{3qk - 2Z) 
Z{2qk - Z){qk - Z) + m^{2qk^ + iZ qk 

m^{3Z-qk) -m^ 



3Z^ 



J", 



(i);8 



ap\Va)k"p{p^, 

Va = {k~pi+ P2)„ {Z^ -2Z + m^) + 
2qk[Z — m?')pi „ 
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4(2) 5 



J(2) 6 




i(2) 7 



7a(l - 75) (P2 - <?2 + "ii) 7p [-Pi - k + mi j 7^ ^ ^ 
7a (1 - 75) {P2 - <?2 + mi)7^ (gi - Pi + wi) 7p ^fa^^ 

((P2 - 92)^ - "if) ((qi - pi)2 - to2-) 
7a(l - 75) (-92 +P2 + ?™i)7A 



Fig. 10. Diagrams for calculation of the process ((Sj and the 
tensor 



i(2) t 



{qi - ky {{p2 - q2Y - ml) 
{f't^ - th") d^^ {qi - k) X 

{(2fc- gi)p5/3M - (fc + 9i)/3.9w + (2'?i - fc)^g/3p} , 
7p (fc - 92 + mi ) 7a(l - 75) (91 -Pi+ "^l)7p 



(fc - 92)^ - m'^^] ((91 -Pi)2 - ml) 



77(1); 9 _ ,v , k'^-n'' X 



[-Z{qk - Z)+ m^{qk - 2Z) + m^] + 
[2g/j2 -3Z qk + + m^{3qk - 2Z) + m^] + 
[iqk^ - Z qk- Z^ + m^{qk + 2Z) - m"*] 



Now we have all the ingredients to calculate ampli- 
tudes of processes ([3|)-([5]). At first let us consider the s- 
channel single top production ([3]), which is shown on the 
Fig.[T]3. We have to calculate N^^. 

From upper and lower diagrams of the Fig.[TT]we have 



77. 



2^2 



rfSf (9)4^'(9)-?l°i. (-91, 0, -92, 0)4°l(pi, m„p2, 0), (83) 



(80) 



Z = f = {pi + k)"^ = (g -p2)^, m = mt. 

And the last one is the amplitude squared of the pro- 
cess which is depicted in the lower Fig. [TUl 

^al PP' '^P^ ' "^1 ' ' '"'i ' 92 , "4 , fc, TiT) = 



h = 



6VDi2{Q^,mt±,Pt) 1-x: 



2 



TT (1 - .T^)2 

(3 + C"2 + (3-C'2):r^), 

C'^ = i?i2(g', P0M2(Q, mt, 0), 
= mw/Q, 

and 



(84) 
(85) 



(81) 



77. 



2-i-3 



,i=l 



Here = ^ = 3 and Cf = {N"^ - l)/(2iV) = 4/3 are 
sructure constants of the group SU(3), color indices are 
contracted with Saa'Sbf in (|8ip . Expressions for Feinman 
diagrams looks as follows 

-^i'j; = v{q2,m'^)A^^lMpi,mi), 



<(9)C'^(9)-^rj'(-9i,0,-92,0)x 
^0^3'{Pi,mt,P2,O,k,K) + 



dwi^- k)d'^^ {q- k)J'^^i^,{pi,mt,p2,0) x 



r(0) , 



J(2) 1 
■^app 



J(2) 2 



i(2) 3 



7p(-9 + "J2)7a(l-75) 


-pi - k + mi^ 


7p 




_ ^^2^ 


(A'2 „ 2pik) 



ft'' 



(82) J-^l^,(-9i,0,-92,0,fc,A-) 
correspondingly, where 
dwi<l)=dwi<l)/i<}^ -mw), d'^iq) 



(86) 



, 9 9 



J(2)4 
^app 



7p (-9 + "^2) 7p (fc - 9 + "^2 j 7a(l - 75) 

(g2 _ ((^ _ fc)2 _ ^^2^) 

7A (fc - 9 + "^2) 7a(l - 75) 

(gi - fc)2 ((q - fc)2 _ ,^^2^ 
(ft'' - t'^e) d^f^ {qi - k) X 

{(2fc- 9i)p5/3p - (fc + 9i)^9pp + (291 - k)^ gpp^ 
7p (fc - 92 + ) 7p ( fc - ^ + 771^ ) 7a(l - 75) 




(92 - A:)2 - TTi!,^ ) ( ((? - fc)2 - m(. 



^ i , Fig. 11. Diagrams for the calculation of the process Q. 
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Fig. 12. Diagrams for the calculation of the process ((4}. 



and for all the calculations wc put = since 

mb/mt < 1. 

For the process ([4]) and calculation of N^^, we have the 
following functions (see diagrams on the Fig. [T^ 



h = -{ -2^/ Di2{Q'^, mt ±,niw ±) x 



(1 - + 2xlu^ (1 + -u^l- + 



xtu' (1 - u')' (I ~ C'^) 



(3 - 2u^) (1 - C"2) + 

3 + 3) - (C"2 + 7) - 2u^ (1 - C- 
xt ( 25 + 3C'2 + 



(23C"2 + 33) - lOu^ (1 - C"2) 

3x1(1 + 2.^) (l-.^)^(l-C'^) 
4(1 + 2u'){l + 2x1(1 - u') + 2xt{l - uY) 



772^2 = (91 -~Pi)d'^^ {qi -_pi) x 



r(0) 



(P2 , 0, -92 , 0)4°l (pi , mt , , 0) , (87) 



^ 288v/7?i2(Q2,mti,Pt) 

J2 = 



1 



{l-xl + 2xl)^-C'^l-xir' 



where C" is the same as in the previous process. 



dwi<l2 ~P2)dw^ (92 -P2) X 



4a' b2, 0, -92, 0) J-^'.l(pi, mt, -91, 0, fc, A-) + 



J-^5,(pi,mt,-9i,0)J-('j,(p2,0,-92,0,fc,AO j>, (89) 



Calculations for the process ([5]) and N^^^ looks as fol- 
lows (sec diagrams on the Fig. [TU)) 



772^2 = (P2) 



(91, 



4a' pp'(Pi."^*> -92, 0,-91,0), (90) 



In 



1 + x^,(l - u^) + ^/Dl2{Q^,mt±,mw ±) 
1 + a;^(l - u2) - A/£'i2(Q^,mt_L,mH--_L) 



., (91) 



where 



C'^ = Di2iQ^, mt ±,mw ±)/£'i2(Q^, mt,mw), 
u = mw /mt, (92) 



77. 



2->3 - 



iP2) 



d"" iqi,n) 



4a' pp, (Pi , , 91 , 0, 92 , 0, fc, 7^ ) . (93) 
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